Coherent transport through a quantum shuttle 
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A shuttle site that crosses through a band of network levels induces currents. If the process is 
re-cycled this can be regarded as "pumping" or "stirring". We show that the analysis reduces to 
that of calculating time dependent probabilities, as in the stochastic formulation, but with splitting 
(branching) ratios that are not bounded within [0, 1]. Our approach allows to address the adiabatic 
regime, as well as the Slow and Fast non-adiabatic regimes, on equal footing. We emphasize as- 
pects that go beyond the familiar picture of sequential Landau-Zener crossings, taking into account 
Wigner-type multilevel mixing due to the crossing shuttle. 



Transport in quantum networks is a theme that 
emerges in diverse contexts, including quantum Hall ef- 
fect Gj , Josephson arrays .quantum computation mod- 
els [3(, quantum internet [J], and even in connection 
with photosynthesis For some specific models there 
are calculations of the induced currents in the adiabatic 
regime for both open and closed systems, so called 

"quantum pumping" ll|-18| and "quantum stirring" [lil - 
|23I | respectively. In the latter context most publications 
focus on 2-level [3, and 3-level dynamics, while the 
larger perspective is rather abstract, notably the "Dirac 
monopolcs picture" B El El HI- This should be con- 
trasted with the analysis of stochastic stirring where the 
theory is quite mature [26M29j | . 

Considering (e.g.) the unidirectional rotation of a 
molecular rotor [27| , it is possibly allowed to be satisfied 
with a stochastic picture [26| that relates the currents, via 
a "decomposition formula" , to rates of change of occupa- 
tion probabilities. Once we turn (e.g.) to the analysis of 
pericyclic reactions [3(| this is no longer possible. In the 
latter case the method of calculating electronic quantum 
fluxes had assumed that the latter can be deduced from 
the continuity equation. Such procedure is obviously not 
applicable for (say) a ring-shaped molecule: due to the 
multiple path geometry there is no obvious relation be- 
tween currents and time variation of probabilities. 

In this work we would like to analyze the following 
prototype problem. Consider a network that consists of 
N interconnected sites, with on-site energies Si, and cou- 
plings Cij. Additionally there is a shuttle (i — 0), where 
the on-site energy £q = u(t) is varied monotonically from 
u = -co to u — oo. Accordingly the Hamiltoinan is 
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Our interest is in the induced current I(t) that flows 
through a tagged connecting bond C a . In the adiabatic 
limit this current is determined by the so-called geomet- 
ric conductance G as follows: 
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FIG. 1: Upper panel: an initially loaded shuttle is crossing 
a band that contains N = 500 levels. The occupation prob- 
abilities of the shuttle (p) and of the levels (q n ) are imaged 
as a function of time. The vertical axis is the energy, and 
the horizontal axis is u(t). We assume star geometry with 
level spacing A=l, and identical couplings c n =3. The dashed 
line illustrates the energy of the lowest adiabatic level. Lower 
panel: plot of q n vs n in several cases. Green line with mark- 
ers - the adiabatic scenario of the upper panel at u — 60. 
Cyan solid line - after decay from a standing shuttle (tt=0). 
Red dashed line - after decay from a fast shuttle (u=5000). 



Here <fi is a test flux through the bond of interest, namely 
C a i-> C a e^, and ^> is the wave- function of the adia- 
batic eigenstate. Specific scenarios to be considered are: 
(i) The shuttle is initially filled with a particle that later 
is transferred to the network; (ii) One of the levels of the 
network is initially filled, and later a current is induced 
via the crossing shuttle. The occupation dynamics in the 
first scenario is illustrated in Fig[IJ which will be further 
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discussed later. For a non-interacting many-body occu- 
pation, results can be obtained by simple summation. 

We shall consider first a Star geometry (inset in Fig[T|) , 
and later a general network geometry. As an appetizer we 
refer to a 3 site model (inset in Fig© . We shall see that 
the integrated current through the C a bond (indicated in 
the figure) is Q = \C a \ 2 /(\C a \ 2 + \C b \ 2 ) for C = 0, while 
Q = CajiCa — Cb) for any Co 7^ 0. It is possible to get 
\Q\ > 1 in the latter case, indicating that a circulation 
current is induced in the system. This result has been 
derived in the past as an approximation, using either the 
"Dirac monopoles picture" , or a two-level calculation. In 
the analysis below the C a are allowed to be large, such 
that numerous network levels are mixed during the shut- 
tle crossing. We also address the non-adiabatic scenario. 

Outline.— We find explicit expression for the cur- 
rent / that is induced either in adiabatic or non-adiabatic 
shuttling process. We first consider "star geometry" (in- 
set in Fig© and later a general network. For the anal- 
ysis we introduce the splitting ratio phenomenology, and 
demonstrate it for a dot-wire "ring geometry" (inset in 
Fig .[3]). We find it essential to distinguish between two 
types of processes: adiabatic crossing and adiabatic meta- 
morphosis. This distinction is important for the under- 
standing of non-adiabatic effects. 

Star geometry.— Let us consider the special case 
of having sites with energies £ n = e n and connections 
Cno = c n , while all the other couplings are zero. An 
adiabatic eigenstate of the system is represented by [aj: 



Using the notation 
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g(E;a,...c N ) = 
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the adiabatic energy is the solution of the equation 
g(E) = E — u. As u is swept from —00 to +00, the energy 
increases monotonically from e n to e„+i, where n is the 
starting level. The normalization factor p = [1 — g'(E)}^ 1 
is the probability to find the particle in the shuttle. For 
the following derivation note that 1/p is a quadratic form 
in c n . Using Eq.© we get after differentiation by parts 
that the current through c n is 
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where q„ = \^> n \ 2 are identifies as the probabilities to 
find the particle in n = 1...N. 

Multiple path geometry.— Needless to say that we 
do not really need Eq. © in order to get the expression 
for G in the case of a star graph. We could simply deduce 



Eq.© from conservation of probability, i.e. from the con- 
tinuity equation I = q n . This is no longer the case if we 
have a multiple path geometry: probability conservation 
alone cannot tell us how the current is split between the 
different paths. Furthermore, we would like to go beyond 
the adiabatic transport formalism, and obtain a formula 
that applies also in non-adibatic circumstances. 

Splitting ratio approach. Let us first see what is 
the expression for G in the case of a general network. It 
is natural to switch from the £i basis to an e n basis that 
diagonalizc the network Hamiltonian in the absence of 
the shuttle. Consequently getting a star geometry with 
c n = ^2b(n\b)Cb- Our interest is in the current through 
a tagged bond C a . We define the "splitting ratio" of the 
current that flows in the nth levels as 



X n [splitting] = 



(n\a)C a 



(n\a)Cg 
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A straightforward generalization of the derivation that 
leads to Eq.© implies that the current through C a is 
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di 
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The physical simplicity of Eq. © suggests that it 
can be derived without assuming adiabaticity, where 
Qn(t) = l^n^)! 2 is obtained from the solution of the 
time dependent Schrodinger equation. Indeed this is the 
case. We just have to remember that quite generally 
I = C o Im[* o (t)** (t)], substitute * a = J2n( a \ n )^n, 
use the definition Eq. © of the splitting ratio, and the 
identification q n = c, i lm['0*'0o]- 

The integrated current.— As the simplest exam- 
ple for the application of the splitting ratio approach 
we consider a process in which a particle has been 
transferred from the shuttle to the wire in the net- 
work that is illustrated in the inset of Fig.[3l In this 
model the splitting ratio of the even-parity levels is 
A n = A + = C a / (C a + Cb), while for the odd-parity levels 
we have A„ = A_ = C a j (C a — Cb)- From Eq.© it follows 
that the integrated current is Q = average(A„), where the 
weighted average is determined by the final occupation of 
wire levels. For an adibatic process, in which the particle 
ends up at the lower wire level, we get Q = Ai. Unlike the 
case of a stochastic transition this value is not bounded 
within [0, 1]. rather it may have any value, depending on 
the relative sign of the amplitudes C a and Cb- However, if 
the process is not adiabatic, the probability is distributed 
over both the odd and the even levels with probabilities 
that are proportional to \C a ± Cb\ 2 respectively. Then 
we get from the weighted average a stochastic- like result, 
namely Q= |C a | 2 /(|C a | 2 + |C h | 2 ). 

It is also instructive to consider the process in 
which a particle that is prepared (say) in an even 
wire-level, is adiabatically transferred, due to shut- 
tling, into the adjacent odd wire-level. Here the in- 
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tcgration over Eq. implies Q = A_ — A+, leading to 
Q = 2C a C b /{\C a \ 2 - \C b \ 2 ). 

The parametric variation of the current. The 

results for the integrated current give the impression that 
the size of the coupling c„ compared with the levels spac- 
ing A is of no importance. But this is a wrong impres- 
sion. Once we get deeper into the analysis it becomes 
clear that the familiar two level approximation for the 
adiabatic current J, requires the coupling c n to be very 
small compared with the level spacing A. Our interest 
below is focused in the case of having a quasi-continuum, 
meaning that the c n are larger than A, hence many levels 
are mixed as the shuttle goes through. 

Before discussing the quasi-continuum case it is useful 
to discuss what happens with the N = 2 wire system if 
the c n are not smaller compared with A. 

Adiabatic metamorphosis.— The heuristic under- 
standing of adiabatic transport is commonly based on the 
analysis of "adiabatic crossing" . Below we shall see that 
we have to take into account an additional prototype pro- 
cess that we call "adiabatic metamorphosis" . The mini- 
mal model for its demonstration requires 3 levels. 

Consider the dot-wire geometry with N — 2. Assume 
that A <C C a ,Cb- Clearly during the adiabatic crossing 
A has no effect, so we can neglect it. The dynamics 
is very simple: there is a dark state \D) = Cb\l) — C a \2) 
that does not participate in the dynamics. The crossing is 
exclusively to \C) = C a \l) + C&|2). Using the same pro- 
cedure as before, one concludes that the splitting ratio is 
Q = \C a \ 2 /{\C a \ 2 + \C b \ 2 ). This looks like a contradiction 
to our previous analysis that predicts Q = C a /(C a — C'b) 
for any cq ^ 0. What is wrong? 

After some thought one realizes that there is a 
second distinct stage in the dynamics during which 
the ground state adiabatically transforms from |C) to 
|e_) = |1) — 1 2). The reason for this metamorphosis is as 
follows: the sites |1) and |2) are directly coupled via Co, 
but also virtually coupled via C a Cb/u. This second-order 
coupling via the shuttle is dominating during the adia- 
batic crossing, but much later, when u becomes larger 
than 7 = |C a C(,/A|, the first-order coupling takes over. 

At first glance one may argue that the metamorphosis 
at such later stage is of no interest: all the probability 
transfer happens during the adiabatic crossing; after that 
the particle is confined to the wire; and can merely "re- 
arrange" itself there. However, this reasoning is mislead- 
ing: we have here multiple path geometry, and therefore 
we can have flow through the shuttle without any transfer 
of probability. In Fig[2]we demonstrate this observation. 
If the metamorphosis stage is taken into account, there 
is no longer contradiction with the general formula. 

Adiabatic mixing.— We turn to treat the dot-wire 
system for JV > 1. For simplicity we focus on some 
energy range where the levels form a comb-like spectrum 




FIG. 2: Flow of the current from the shuttle to an N = 2 
wire through the C a bond. The parameters are C a = 19, and 
C'b — 15, and Co = 1. The thickest line is the exact adiabatic 
result for G(u). The thinner and the thinnest lines are I/ii 
for ii = 2 and for u = 50. The left and right vertical lines in- 
dicate the shuttle- wire crossing point, and the metamorphosis 
point, with separation 7 = 285. During the adiabatic meta- 
morphosis a current is flowing through the distant shuttle. 



with level spacing A, hence from the definition Eq.((4]) 
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where c± = prefactor x (C a ± C&) are the couplings to 
the even and odd levels respectively. For simplicity of 
presentation we absorb the prefactor into the definition 
of C a and Cb Q • Of interest is the case of having a quasi 
continuum, meaning that the couplings c„ are larger com- 
pared with A, hence a two level approximation is out of 
the question. We shall see that the role of 7 in the N = 2 
analysis is taken by 7 = sinf T, where 



u + 



sin(0) = 



(9) 



(here and below we assume starting with even-parity 
level). With Eq. ((8|) the secular equation g(E) = E — u 
becomes a quadratic equation for tan(), and can be solved 
explicitly. Then it is possible to get an expression for the 
shuttle occupation probability: 

p(u) = [l-g'(E)}- 1 = A-L[u-E;T,e] (10) 

where the distorted Lorentzian L [x; F; 9] is 
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In the expression above E is the energy in which the 
particle has been prepared, and it can be regarded as a 
constant. Some further straightforward algebra leads to 



G{u) 



au 



(n\a) 



d_ 

du 



C„ 
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c 2 sm 2 (ip(u)) + c 2 _ cos 2 (<p(u)) 
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FIG. 3: The variation of the occupation probabilities as a 
function of u for an adiabatic shuttling process. The particle 
initially has been placed at n = 250. The couplings are C a — 6 
and Cb — 4. The red thick line is p. The other solid lines are 
Q250, and 9251- The dashed lines from up to down are 5249 
and ?253 and q 2 47 and 17252 • 
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FIG. 4: Flow of the current from the shuttle to the wire 
through the C a bond, for the same scenario as in Figj3] The 
parameters are indicted in the legend. In the left panel A = 
200, hence a two level approximation is satisfactory. In the 
right panel A = 1, and consequently the line shape reflects a 
multi-level crossing. 



where ip(u) ~ (tt/2)(E/ A) has a ir/2 variation as u is 
raised from — oo to +00, hence Eq.(|12j) is a step-like func- 
tion that varies from A + to A_ as had been anticipated 
by inspection of Eq.([7]). 

The functions p(u) and G(u) are plotted in Figj3] and 
in FigU] In the latter we contrast with the c„< A case, 
for which the dynamics can be regarded as a sequence of 
two- level crossings. 

Regimes.— The results for the integrated current give 
another wrong impression: it looks as if we are dealing 
with two regimes: either the process is adiabatic or non- 
adiabatic. A more careful inspection reveals that de- 
pending on u we have 3 regimes: Adiabatic, Slow and 
Fast. For star geometry with comb-like quasi continuum 
of levels, the Slow regime is defined by the condition 



< u < r 2 , 



T = 2tt- 



A 



(13) 



For simplicity we assume here comb-like quasi contin- 
uum with identical couplings c„ = c. The left inequality 
in Eq. (|13p means that the adiabatic condition is violated, 



while the right inequality implies that a first-order per- 
turbativc approximation is violated as well. The iden- 
tification of this intermediate Slow regime parallels the 
notion of Wigner or FGR or Kubo regime in past studies 
of time dependent dynamics 

Some illustrations for energy spreading are presented 
in FigJTJ If c < A the transport of probability from the 
shuttle to the wire would be described using a two level 
approximation. But the illustration in the upper panel 
assumes c > A, hence many levels are mixed within a 
parametric range T. The time during which this mixing 
takes place is T/ii. In the opposite limit of Fast shuttling, 
which we further discuss below, the decay time of the 
probability to the quasi-continuum is 1/F. 

Non adiabatic spreading.— The calculation of / 
in the non-adiabatic regime requires knowledge of q n (t) . 
For star geometry this calculation is a variant of the 
Wigner decay problem, and hence can be solved ana- 
lytically: instead of a fixed level that decays into a quasi- 
continuum we have a moving shuttle. The usual text- 
book procedure is followed [31| leading to the equations 
<9 t 4?o = [-m(t)-(r/2)]* , andt9 t *„ = -ie„4?„-ic„4?o- 
With u(t) — ut one obtains the solution 



Qn(t) 



i / dr cxp 





ie n r 



(14) 



By inspection one observes that going from the Slow 
to the Fast regime, the spreading line shape changes 
from Lorcntzian-type to Fresnel-type, as illustrated in 
the lower panel of FigQ] 

Discussion.— We have found, using elementary con- 
siderations, without the need to rely on a complicated 
transport formalism, that it is possible to replace Eq.([2]) 
by the general expression Eq.©, that holds both in adia- 
batic and non-adiabatic circumstances. Hence the prob- 
lem of calculating currents is reduced to that of calculat- 
ing time dependent probabilities q n (t) as in the stochas- 
tic formulation [2f3 |. It is important to realize that the 
"splitting ratio" Eq.© unlike the stochastic "partition- 
ing ratio" is not bounded within [0, 1]. This observation 
has implications on the calculation of "counting statis- 
tics" and "shot noise" [HHIl. 

In the analysis we have emphasized aspects that go be- 
yond the familiar two-level approximation phenomenol- 
ogy. Whenever we have quasi-continuum the network 
levels are scrambled during the shuttling process: this 
is what we called "metamorphosis" or "mixing". This 
scrambling is reflected in the time dependence of the in- 
duced currents. 

Finally, we realize that the non-adiabatic dynamics is 
characterized by a stochastic-like result for the integrated 
current Q. But the detailed temporal variation of the 
current / has different features depending on whether 
the shuttling process is Slow or Fast. 
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Eigenstates for star geometry.— For a star geometry we have 



H = 



Cl £i 

C2 e 2 



-N 







(15) 



\c N ... e N j 
The eigenstates satisfy the following set of equations: 



i*0 + X] C n*n = 



From Eq. ((T7)) we get 



Cn^o + e„^ n = E$ n , n=l,2,...,N 



(16) 
(17) 

(18) 



Substitution to Eq. (fTB"|) gives the secular equation g(E) = E — u, as stated after Eq. Q. From the normalization 
condition we get the expression for p, namely 



N 



ii-g'(E)]- 1 



(19) 



Calculations of p for a dot-wire geometry.— The energy levels of the wire that has length L = N + 1, are 
e n = — 2cq cos(fc rl ), where k n = {-k/L)ti. The respective couplings to the shuttle are 



L 



1/2 



sm(k n ) 



(C a ± C b ) 



(20) 



where the ± reflects the parity of the level. We focus on levels with energy e„ ~ E, such that their spacing A can be 
regarded as constant. Then it is convenient to absorb the factor [(2/N) 1 / 2 sin(fc)] into the definition of C a and Cb- 
The couplings c n = c± to the energy levels in the sum Eq.Q are distinguished by their odd/even parity, hence after 
summation we get two terms: 



2 ,1 E\ 2 / E 



(21) 



The secular equation g(E) = E — u becomes a quadratic equation for tan(), and can be solved explicitly. 



COt 7T 
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(22) 



where the ± refers to the parity that is alternating for subsequent levels. Then it is possible to get an explicit 
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expression for the shuttle occupation probability 
P = [l-g'(E)}- 1 = 
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where the distorted Lorentzian is defined as follows 
L [x; T; 6] = 



1 



smtf x 



rr 2 +cos 2 <9 (r/2) 2 ] 1/2 



cos 2 6 (r/2) 
7T J x 2 + cos 2 6> (r/2) 5 



Note that for = (c + = c_) and for = 7r/2 (c_ = 0) it becomes the standard Wigner Lorentzian. 
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Calculations of G for a dot-wire geometry.— The calculations of G involves similar sums as in the p calculation. 



leading to Eq.(fT2 



G = 



du 
C a - 



= C a 
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-E 
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\ 7T ^ + sin^^) 
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Calculations of G for a 3 site geometry.— Here we show how to to re-derive the exact result for G(u) in the 
case of a 3-site system [an using the splitting ratio approach. For such system the 2 wire levels e± = ±cq have 
coupling c± = (C a ± Cb)/V2 to the shuttle, and accordingly A± = C a /(C a ± C&). Using Eq.Q we get 



G = — 
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